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The propagation and non-linear interactions of magnetohydrodynamic waves are considered in 
the force-free limit, where the inertia of the conducting matter which enforces the MHD condition 
E ■ B = can be neglected in comparison with the inertia of the electromagnetic field. By extending 
the analysis beyond the WKB approximation, we are able to study the non-linearities induced by a 
gravitational field. We treat the perturbed electromagnetic field as a fluid of infinite conductivity. 
, We calculate the scattering of a torsional (Alfven) wave by a gravitational potential, and demonstrate 

■ a nonlinear coupling with a compressive (fast) wave which is second order in the amplitude of the 

Alfven wave. In a cylindrically symmetric spacetime with slow rotation, the coupling is second 
£Sj . order in g t( p and first order in the amplitude of the wave. We also give a fresh analysis of the 

I ■ non-linear interactions between compressive and torsional waves in Minkowski space, with a focus 

on the relative strengths of their three- and four-mode interactions. In contrast with non-relativistic 
magnctofluids, the effects of compression are always present. In the case of colliding fast waves, 
£ — ■ a net displacement of the field lines across (at least) one of the colliding wavepackets is shown to 

H ' have a strong effect on the outgoing waveform, and to have a qualitatively different interpretation 

than was previously suggested for colliding Alfven waves. Finally, we show how spacetime curvature 
7-H ■ modifies the collision between two torsional waves, in both the weak- and strong-field regimes. 
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I. INTRODUCTION 

Oh| Magnetohydrodynamics can be studied in the ultra-relativistic limit, where the energy density in the conducting 
matter which enforces E • B = is much less than in the electromagnetic field itself. The Euler equation reduces 
to the force-free equation J^F^ = 0, which provides a useful starting point for investigating some aspects of the 
dynamics of accreting and outflowing matter around compact stars, and classical gamma-ray bursts (e.g. Q], 0, 0, 
0). In this regime a uniform magnetofluid supports two modes: incompressible torsional waves analogous to Alfven 
waves in non-relativistic magnetohydrodynamics (MHD) , and a compressible wave analogous to the fast mode. Both 
of these modes have exact non-linear solutions in a uniform magnetic field in Minkowski space, which means in turn 
that there is no spontaneous decay of one type of mode into the other. 

This fundamental property of a relativistic magnetofluid changes in the presence of a strong gravitational field. We 
focus in this paper on the effect of spacetime curvature on the interactions between torsional and compressive MHD 
modes. There is a non-linear interaction if the background spacetime is static, which appears at linear order if it is 
slowly rotating. We analyze the effect first in the weak-field regime, and consider both spherically and cylindrically 
symmetric gravitational fields. As a particular example where strong-field effects are important, we also consider the 
case of a uniform magnetofluid immersed in the spacetime of a black string. 

Our main goal here is to explore the fundamentals of the interactions between MHD waves and gravitational fields. 
Although the scattering of gravitational waves, scalar waves, and vaccum electromagnetic waves have long been studied 
in black hole spacetimcs 5] , analogous work on force- free magnetohydrodynamics has focused on the short- wavelength 
limit 2]. There are two well-established astrophysical contexts in which the dynamics may be well described, in a first 
approximation, by a relativistic force-free fluid: the magnetosphere of a bursting magnetar (e.g. 0); and the polar 
regions of an accreting black hole [l|. We restrict ourselves, in this paper, to the simplest case where the background 
magnetofluid does not maintain a large-scale current (as it does surrounding a Kerr black hole). 
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In a dynamic situation, a fluid formulation of force-free magnetohydrodynamics provides a physically transparent 
set of variables. Two such formalisms have been developed, by Achterberg Q and Thompson & Blaes [3|. These are, 
in fact, dual versions of each other, and are connected by interchanging the dynamical role of the Lagrangian fluid 
coordinates and the Eulerian coordinates of the background spacetime. In reviewing these formalisms, we generalize 
the second to an arbitrary curvilinear coordinate system. 

Another of our goals is to extend the discussion of non- linear mode interactions in force- free MHD, begun in 
We examine in more detail how obliquely propagating Alfven modes and fast modes will perturb each other, including 
the effects of spacetime curvature. The effect of field line wandering on mode collisions has received considerable 
attention in recent years, in the case of non-relativistic and incompressible Alfven turbulence. We examine its effect 
on the collisions between compressive (fast) modes in a relativistic magnetofluid in flat Minkowski space. In contrast 
with the case of colliding Alfven waves, the effect of a net field-line displacement cannot be expressed in kinematic 
terms through a zero-frequency component of one of the colliding waves. 

The sections are organized as follows. In Section II, we given an overview the electric and magnetic lagrangian 
fluid formulations of force-free MHD. As a toy example, we construct the solution for the uniform magnetic field in a 
spherically symmetric spacetime (with non- vanishing Ricci tensor). In Section III we review the normal modes of a 
uniform, force-free magnetofluid. Section IV is devoted to the scattering of a torsional wave in a curved spacetime: 
we first consider the propagation of such a wave in a general, cylindrically symmetric spacetime; and then consider 
the scattering of such a wave in a shallow, spherical gravitational potential. Section V generalizes the results of the 
previous section to the case of two colliding, axisymmetric torsional modes in a gravitational field. The final Section 
VI considers the more general case of collisions between non-axisymmetric MHD modes in Minkowski space, focusing 
on the three-mode couplings of two colliding Alfven and fast waves. 



II. LAGRANGIAN FORMULATIONS OF ULTRARELATIVISTIC MHD 



The equations describing the behavior of ideal fluids coupled to electromagnetic fields in curved spacetimes are 
inherently very complicated. A significant simplification is obtained by assuming that the spacetime metric is fixed; 
i.e., that the mass-energy of the fluid can be neglected in comparison with that of the source of gravity. In this paper 
we make a second simplification: the role of the charged matter fields is restricted to a source of electric currents 
which cancel off the electric field in the fluid rest frame. The mass-energy of these matter fields is, in other words, 
negligible in comparison with that of the electromagnetic field. 

We begin by summarizing the two Lagrangian fluid formulations of relativistic, force-free MHD. 



A. MHD equations in the force-free limit 



The covariant Maxwell Equations are 



and 



V F + V F 4-VF — 8 F 4- 8 F 4- 8 F =0 



(1) 



(2) 



When the only force acting on the charged matter in a local incrtial frame is the Lorcntz force, the dynamics of the 
electromagnetic field is described by the force-free equation 



in combination with 



Here 



T^F ——V F^ a F — f) 

47T 



F^F^ = 0. 



a/3 



(3) 
(4) 
(5) 



is the dual of the electromagnetic field tensor, expressed in terms of the antisymmetric symbol with unit norm 
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B. Electric lagrangian variables 

Achterberg Q developed a variational principle for relativistic magnetohydrodynamics in which the Lagrangian 
coordinates Xq of the fluid are the field variables, moving in a fixed Eulerian space x^ . In the ideal MHD limit, where 
the electric field vanishes in the rest frame of the fluid, the Lie derivative of the Faraday two-form vanishes, £ U F = 0, 
along the direction u of fluid motion. One thereby obtains a simple relation between the background (reference) state 
and the dynamic state of the magnetofluid. Solving for the Xq allows one to obtain all relevant information about the 
system. 

The Faraday two-form is related to its background counterpart by 

_ dXQ a dXpf* (Q) 

~ dx» dx» ' W 

In the force- free limit, the action is 

S = - J ^x-^—gg^g^F^F pa . (7) 
By extremizing this action, the equation of motion 

J^F^ = (8) 
can be obtained. Written in terms of the coordinate fields this is 

dx Oj?(o) 8 ( r^^^r9x 7 Q dx s {a) 



dx» av dxPy 99 9 dx°dxT^)- u - W 

Here, we have eliminated a factor of Oxq/Ox 1 ' from F^. The Bianchi identities are automatically satisfied as long 
as they are satisfied by the background field. Once the equations for these coordinate fields have been solved, the 
electromagnetic field is obtained from Eq. ©. 



C. Magnetic lagrangian variables 

Two new MHD formalisms in the extreme relativistic limit were recently developed by Thompson and Blaes 
The first of these is related by a duality transformation to the action presented by Achterberg, and is of interest here. 
In this formalism, the dual tensor F (rather than F) plays a central role and the dynamical fields are the perturbed 
coordinates x M (rather than the Xq). In this section, we generalize the treatment of Q to include gravity. 

Once the fields x^ are solved, the electromagnetic fields can be found from the transformation 



pnu = V-9(xo) 9xf dx v ~ ap 



(10) 



Here 



(11) 



is the Jacobian of the transformation x$ — > a; M from Lagrangian to Eulerian variables, and the factor of y/go/g is 
necessary in a curvilinear coordinate system. We show in Appendix lAl that the force-free equation J^F^ — may 
be derived from the action 



S' = 



J d 4 xoL' = J d 4 x—' 



-99w9u<x 



d Xq 



-gg^ P gua 



(12) 



Along the way, we obtain the equation of motion 



F d 



a 

dx £ 



gpsg-wF' 



Pi 



l 

2 F 



5e 



dx^ 



gpsg^eF 



0i 



(13) 
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In cases where the background magnetic field asymptotes to a constant at large distances, we can use the simple 
background 

= BoVffi - Sffi); = 1, (14) 



— "0V U 0"3 "3 

and define proper time t = to and distance a = zq. The equation of motion Ijl3|l then takes the form 

daf_ fdx^d_ _ dx^d_\ f dx^dxf\ _d_^ _ 

9xq \ <9r 9cr da dr J Se \ dr da J dx% &e 



where 

9iia9vfi ( dx a dx@ dx 13 dx a 



- ^=t ( -ir-^r - -iz-^z ) = -5V (16) 



gJ4 \ dr da dr da 

Although Eq. I|13ll has four components, it is straightforward to check that the longitudinal components fj, = 0, 3 
vanish. Physically, this is equivalent to the statement that MHD waves have two independent transverse polarizations. 
Notice also that it is possible to simplify Eq. (|13|) further by choosing a time slicing in which 

V^gJi = 1. (17) 

The cost of this gauge choice is that one can no longer identify the background time coordinate t with the proper 
time r of the magnetofluid. 

D. Uniform magnetic field: Spherically symmetric spacetime 

There is a simple prescription for constructing a uniform magnetic field ( — 0) in a vacuum spacetime with an 
axial killing vector £ M = g^^. One identifies Q 

4, = \bo£* = \bq94>^ (18) 

Thus Fuv takes exactly the same form outside a black hole of vanishing charge and spin, in the standard Schwarzschild 
coordinates (g^ — r 2 sin 2 6), as it does in Minkowski space: 

F r4> = B r sin 2 0; F e<j> = B r 2 sin cos 0. (19) 

More generally this method is of limited utility, since it depends on the vanishing of the Ricci tensor of the 
background spacetime. To give a further illustration of the lagrangian fluid method, we show how it may be used to 
construct a uniform magnetic field in a general spherical spacetime. The metric may be written as 

ds 2 = g tt dt 2 + g rr dr 2 + r 2 (d9 2 + sin 2 0d(f> 2 ). (20) 

We imagine that the currents supporting the magnetic field are flowing far outside the self-gravitating mass, which 
itself is not electrically conducting. Then 

Ffiv = (d^RodM - d^ d u Ro)F^ = i(0„.Rgft,0 - d^ d v R 2 )B . (21) 

Here i?o = r sin O is the cylindrical radius and FjL = B R, is the field strength that would be supported by the 
same currents in flat space. The vanishing of 

47r^J = d r (V=gg^g rr F^ r ) + d e (V=gg^g 90 F^) (22) 

combined with the separation of variables 

rl sin 2 0q = K(r) T{0) (23) 

gives the coupled ordinary differential equations 

d ( I ch7dlZ\ sui0 d ( 1 e?T\ 

' ' 1 (24) 



7Zy '— gttgrr dr \y g rr dr J T d6 \sm0 d9 
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Because the equation for T is the same as in fiat space, we may write 

sin0o = sni0 (25) 



and 



where 1Z is a solution to 



d f l—gtt d1Z 



^dr^^j-'^^n^O. (27) 

If the Ricci tensor of the background spacetime vanishes at large radius, then the appropriate solution to l|27|l must 
match onto lZ(r) ~ r 2 at large r (where —gtt = 9rr = 1 — 2GM/r). 



E. Uniform electromagnetic field: Cylindrically symmetric spacetime 

We will be concerned with the propagation of MHD waves in cylindrically symmetric geometries of the following 
form 

ds = gttdt + g rr dr + g^dcj) + g zz dz + 2g tc / ) dtd(j). (28) 

When g^ ^ 0, the frame dragging of the background z-magnetic field generates a radial electric field. Thus, the 
background electromagnetic field is a stationary, axisymmetric, and z-independent solution to the equation = 
with non- vanishing components and , 

d r (V~gg tt g rr F^ + V~gg^g rr F^) = o ; (29) 



d r (V—gg^g rr F$ + V^gg^g rr F^) 

Since the electric field vanishes in the static limit, one has 



= 0. 



(30) 



p (o) = 9p F (o) 

rt tt r<t> 

go 



and 



Bn 



go rrgo 4>4> 

V-go 



B . 



(31) 



(32) 



[Here 50 = s(a;o)-] 

We will investigate perturbations to the fields, subject to the constraint F^F^ — 0, in spacetimes with varying 
degrees of curvature. It is useful first to consider the case of a nearly flat and non-rotating cylindrical spacetime. In 
general, static spacetimes which are symmetric about the z-axis may be written in following form 9] 



ds 2 



„2A 



dt 2 + e 2 ^- x \dr 2 + dz 2 )+r 2 e~ 



2A j i2 



(33) 



The function A satisfies Laplace's equation in cylindrical coordinates, and v is second order in A. In the weak field 
limit this metric reduces to 



and 



ds 2 = -(1 + 2X)dt 2 + (1 - 2A)(dr 2 + dz 2 ) + r 2 (l - 2X)d(/) 2 



-G I dV P{X ' ] 



\x — x'\ 



(34) 



(35) 
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becomes the usual Newtonian potential. The background fields (|32fl and (|31|l now reduce to 



F^=(l-2X)B r; 



= 0. 



(36) 



A cylindrically symmetric spacetime with some attributes of 3 + 1-dimcnsional black holes can be constructed from 
the product of 2 + 1-dimensional black hole (sitting in the r-cj> plane) and an infinite line segment along the z-axis. 
The 2 + 1 black hole solution was originally found by Banados, Teitelboim and Zanelli its properties are reviewed 
by Carlip . The key element in its construction is the introduction of a negative cosmological constant A = — -m , 
so that the spacetime asymptotes to anti-de Sitter space at large radius. The BTZ metric is 



ds 



r 



4r 2 



dif 



r 

- M+ P 



4r 2 



dr 



2^ dt 1 " r 



(37) 



or equivalently 



ds 2 = -(ZR)dt 2 + [ZR)- 1 dr 1 



J 

2^2 



dt 



dz 



(38) 



where 



-M 



T 



(39) 



and 



R=l- 



= 1 



J 



9ttg^4> 



AZr 2 ' 



(40) 



(The variable M, being a mass per unit length, is dimensionless in units where G = c = 1.) The background Maxwell 
fields in this spacetime are 



F. 



(0) 



B r 
ZR' 



(41) 



These expressions reduce to 



p(0) _ B o r p (0) _ n , 4 o\ 

r( t> ~ r*/P-M' rt ~ ( ' 
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in the static limit (J = 0). Minkowski space is recovered by taking M — > — 1 and 1 2 — > oo 



III. FORCE-FREE MHD WAVES IN CURVED SPACETIMES 



A. MHD modes: Force-free limit 



We first review the normal modes of a uniform magnetofluid in flat Minkowski space, before considering their 
interactions and the effects of spacetime curvature. The background magnetic field is B — Bqz, or equivalently 
F ri p — Bar in cylindrical coordinates. The fluid supports two distinct modes in the force-free limit: the fast mode, 
which involves a compressive disturbance and has an isotropic dispersion relation to = |k|; and the Alfven mode, 
which is incompressible and has a group velocity directed along the background magnetic field, duj/dk z = ±1. There 
is no slow mode. 

The simplest way of seeing that the fluid supports only two normal modes is to note that the perturbation to the 
Faraday tensor 



SF — F - 



(44) 
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can be expressed in terms of two variables 5r and 5<f>. To first order in the perturbation, we have 



SF $ = S- (F<$5- 



rip Q r \ M r<f> " ) + f)^ rep 



d5r 



dSr (o) 



z0 Q z r<f> ' 



and 



, (i) _ dS(f> (0)j 



«5ir« = -^F 



dz 



sf£ = 0. 



(45) 



(46) 



As expected, the component of the electric field parallel to Bo vanishes to linear order. The absence of the slow mode 
can be related to the invariance of Ffh) under reparameterizations of the z-coordinate. The full expressions for 8F^ V , 
written in terms of the Lagrangian fluid variables, are collected in Appendix A. 

More generally, the Alfven mode and fast mode involve perturbations to both fluid coordinates. These perturbations 
are related by the condition of incompressibility 



1 d(r8r) d8<f) 



(Alfven) 



r dr 

in the case of the Alfven mode; and the condition of vanishing torsion 

dSr d{r 2 54>) 



dr 



= (Fast) 



(47) 



(48) 



in the case of the fast mode. The fluctuating current associated with the field perturbations l|45|l. I|46|) is 

Av^—gJ^ = d r (V~99 U g rr 5F^) + 9 (V^flV*^) 5 ( 49 ) 



4ttV^9J 2(1) = d r L 



^g zz g rr 5F^ + 9 (y^a"a**SF^ ] 



-gg-g^of^' ; 



(50) 



47^^ 



-gJ r ^ = d t (V=99 rr 9 tt SF$) + d z (^g"^6F$) + d (^gg rr g^5F$) ■ 



and 



47rV^ (1) = d r (y^g^g rr SF^) +d t {y=g^g u 6F^)+d t (y=gg**g"6Fg) . 

Expressed in terms of the coordinate perturbations, and specializing to Minkowski space, this becomes 

AnrJ^ = B d t [d r (r 2 60) - d^Sr] ; 



(51) 



(52) 



(53) 



Airr J z « = -B d z [d r (r 2 S4>) - d^Sr] 



(54) 



and 



47rrJ r 



B r 2 -d, 



1 



B i -d?Sr + d 2 Jr + d r 



-d r (rSr) + d^Stj) 



-d r (rSr) + d$l 



(55) 



(56) 



The dynamical equations for the two normal modes are recovered by imposing the force-free condition. To linear 
order one has 



r(i)p(0) _ _ 70(i) p(0) 



(57) 
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hence J r ^ = J<t>(X) = o. In the case of the Alfven mode, applying the constraint (|47ll gives the one-dimensional wave 
equation 



- d 2 t 5r + d 2 Sr = 0; 
d 2 5<j) + d 2 8<p = 0, 



with a vanishing longitudinal field perturbation 



SF. 



(i) 



0. 



and (in general) a non-vanishing longitudinal current 

jt(D ) ^o. 



(58) 



(59) 



(60) 



A force-free magnetofiuid supports torsional waves of both helicities. It should be noted that, in this limit, there is 
no basic distinction between the two helicities, as there is a non-relativistic cold, magnetized plasma. 

In the case of the fast mode, one instead finds that the current vanishes entirely and the coordinate perturbations 
obey three-dimensional wave equations, 



d 2 S<P + d 2 z 6(f> + \d r [rd r (r 2 <^)] + 



= 0; 



d?Sr + d 2 Sr + d r 



d r (rSr) 



+ ^dpr = -d^8<j). 



By applying the operator r d r r to the first equation, one sees that an axisymmetric field perturbation 



5B Z 



SF, 



(i) 



<f>T 



F, 



(0) 



<f>r 



1 9 (A \ 
-— (rdr) 

r or 



(61) 



(62) 



satisfies the (axisymmetric) cylindrical wave equation. 

Thus, the fast mode is equivalent to a vacuum electromagnetic wave superposed on the uniform background magnetic 
field, but with its polarization restricted to <5E ■ B = 0. This way of constructing the fast mode also makes it 
immediately apparent that infinite plane waves (restricted to this single polarization state) are another general class 
of solutions. 

Inspection of Eqs. I|55(l and l|56(l shows that another projection, 



1 d8r 85(f) 
r 2 d(j) dr 

allows the equations of motion for Sr and S<f> to be entirely separated, 

- 8 2 6<j) + d 2 z 8<j) + \d r (r 3 d r 50) + \dl5(t> = 0; 



(63) 



d?Sr + d 2 Sr + d r 



d r (rSr) 



1 



dlSr = 0. 



At the same time, the longitudinal currents do not vanish, 

2tt 2tt 



(64) 



(65) 



B. Decomposition into normal modes in a curved spacetime 



We have seen that, even in Minkowski space, there is an ambiguity in the definition of the fast mode. Axially 
symmetric perturbations decompose directly into torsional modes, which are supported only by a perturbation Sip, 
and compressive modes, which are supported only by a perturbation Sr. In the more general case of non-axisymmetric 
perturbations, one can still define unambiguously an incompressible mode ()47|l . But one is faced with a choice either 
of defining the fast mode to be current free through Eq. (|48|l , in which case the equations of motion (|61|l for Sr and 
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8<p do not entirely separate; or instead of using the alternative projection (|63|l which separates these equations but 
leaves a non- vanishing longitudinal current (|65[) . 

This ambiguity remains when one considers the case of a magnetofluid in a gravitational field. We focus here on a 
general static {gt<f, = 0) and cylindrically symmetric spacetimc. In this case, the longitudinal currents (|49|l and 15UI) 
cannot both vanish, as is seen by writing them in terms of the coordinate perturbations: 



and 



s v 9 at 



B dz 



d_ 

dr V g tt 
d ( 9u 



dr \g 



g rr dSr 
9tt d(j) 



g rr dSr 



fjz 



(66) 



(67) 



In these expressions, gtt and g zz have different dependences on radius r, and we have made use of the diagonality of 
the metric. We conclude that, to linear order in the field perturbation, non-axisymmetric excitations of the fast mode 
are current-carrying. 

It is, nonetheless, possible to choose a generalization of the projection (|63[1 which does allow the equations of motion 
for dr and 5<f> to decouple at linear order. From Eqs. I|66|l and (|67|l . the obvious choice is 



g rr dSr dScj) 

~5T a - = 0. (Fast) 

g^ d(p or 



Then Eq. lj5l21) becomes 

1 d 2 Sr 



1 d 2 Sr 



F, 



(0) 



1 d 2 Sr (o) 



6> 



,5tt 5t 2 ' g zz dz 2 1 g ^ d(j) 2 

or, upon substituting expression 132(1 . 



_ 5 3 



1 d 2 Sr 1 <9 2 <5r 



1 d 2 Sr 



gtt dt 2 g zz dz 2 
The longitudinal currents then are 



J_d_ 

grr dr 



f^gdr 

g<p<j>grr dr 



g^g rr dr 



(8rF { r , 



(<>) 



0, 



dr 



grrg&tf} 



= 0. 



So 



dr \ g tt J dt 



B n 



gJ z 



dr \ g 



dz 



(68) 



(69) 



(70) 



(71) 



The perturbation 
is defined by 



can be recovered from 5r using the constraint l]68p. Similarly, the incompressible Alfven mode 
V~g d ( grrg^ 



grrg<p<p dr 



°£ = 0, (Alfven) 



(72) 



to linear order. 

For example, the wave equation for a purely compressive disturbance in a BTZ black string spacetime is 



1 d 2 Sr 1 d 2 5r 



1 d 2 8r 



Z 2 dt 2 Z dz 2 r 2 Z 



d 


Z d 


( rSr\ 


dr 


r dr 





= 0. 



(73) 



(As before, Z(r) = r 2 /£ 2 — M .) When d z 5r = d^Sr — 0, this equation can also be re-expressed directly in terms of 
the field perturbation 



SB Z 
~B* 



Z d frSr 
~rdr \~Z~ 



as 



(F_ f SB 2 
dt 2 \~B^ 



Z_d_ 

r dr 



Zr dr { B z 



0. 



(74) 



(75) 



More generally, the field perturbations have some dependence on z, and it is more useful to work with the coordinate 
perturbation. 
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IV. SCATTERING OF TORSIONAL WAVES BY SPACETIME CURVATURE 



We now turn to an important, but subtle, difference between the compressive (fast) and torsional (Alfven) modes 
in a curved spacetime. We focus on a uniform magnetofluid in which the background current vanishes. Both modes 
are exact non-linear solutions to the force-free equation in Minkowski space. This property is retained by a purely 
radial oscillation (6<j) = 0), in any static cylindrically symmetric (or spherically symmetric) spacetime. The resulting 
fast mode is equivalent to an electromagnetic wave polarized in the (^-direction, with <5E • B = 0. Hence it can have 
arbitrary amplitude. 

An axisymmetric torsional wave is not, by contrast, a non-linear solution to J^F^ U = 0, because the wave speed 
varies with radius. We now show that, in a static spacetime, such an Alfven wave will excite a compressive motion 
transverse to the background magnetic field. The amplitude of this compressive motion is second order in the 
amplitude e of the Alfven wave itself. This effect is somewhat analogous to the non-linear coupling between a fast 
mode and an Alfven wave that is propagating through a non-relativistic magnetofluid with a gradient in density and 
in Alfven speed (e.g. |12| . 0|). In a spherical spacetime, we show that a torsional wave originating at infinity is 
scattered into a spherical outgoing compressive wave. 

To this end, we expand the field and currents perturbations in powers of e <C 1: 

Fia , = Fjg) + 6F$ + 5F$ + 0(e 3 ); (76) 



J" = + 5J^ 2) + 0(e 3 ). (77) 

Here, SF^ = 0(e N ), etc. 

When considering the effects of spacetime curvature, we also restrict ourselves to axisymmetric perturbations. In 
a static spacetime, this means that 5r and 8<f> decouple to linear order. Such a linear coupling is re-established in a 
rotating spacetime. 



A. Torsional waves in a cylindrically symmetric spacetime 

1. Static spacetime (g t $ = 0) 

Torsional deformations of the magnetic field are supported by electric currents which, to leading order in e, propagate 
along the background magnetic field. We first derive the analog, in a curved cylindrical spacetime, of the linear wave 
equation l|58|). With this in hand, we find solutions to J^F^ = to second order in e. 

The force-free equation implies J r ^ = to linear order (Eq. |57|1. Combining expression (|51|) with the with the 
Bianchi identity dtSFrl = d z 5F$ gives 

g^(sF^)+g^(6F^)^0. (78) 

The electric field 5F^ satisfies the same wave equation. The solution is 

<5F± = 6F±(z t v,t); 6F± = 6F±(z T v„t). (79) 

where 



(80) 

The two fluctuating field components are related by 

5F± = Tv z 5Ff r . (81) 

The net result is that the Alfven wave propagates along the z-axis, to leading order in e, with a phase speed equal to 
that of a null geodesic that lies tangent to the background magnetic field. 

An Alfven wave front propagating along the background magnetic field will be bent as the result of the dependence 
of gu (and v z ) on radius. This, in turn, forces a radial oscillation of the magnetofluid that is second order in e. Indeed, 
inspection of the force-free equation 

J^SF^ + rMSFjU + J«*>Fg> - 0, (82) 
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shows that the first two terms generally do not cancel. Using Eq. (|81|l to eliminate SF tr in terms of SF zr , one obtains 



rr zz t 

9 9 



■■(SF^) 2 d r (\nv z )+4nJ^F^ = 0. 
The radial fluctuation of the magnetofluid supports the additional current perturbation, 



4ttV 3 5J 0(2 



-Ba 



g rr d 2 Sr g rr d 2 5r d 
gtt dt 2 g zz dz 2 dr 



9 d ( grrg<t>4> 



g<f>4>grr dr 



(83) 



(84) 



Let us first consider the special case of a nearly flat cylindrical spacetime. In such a background, v z — 1 + 2A and 
Eq. (JSSl reduces to 



to lowest order in A. Re-expressing this equation in terms of the perturbations Sr and S<p gives 



d 2 Sr d 
~dt? + ~d~r 



ld_ 

r dr 



(rSr) 



d 2 Sr 
dz 2 



= -2r 2 



dX fdScj) 



dr V d 



(85) 



(86) 



The analogous expressions for the BTZ black string spacetime are not much more complicated. In this case, 
v z = Z 1 / 2 = (r 2 /£ 2 — A/) 1 / 2 , which vanishes at the horizon r# = M x l 2 l. The force-free equation reduces to 



\{5F^fd r Z + 4^ 2 >4° r ) =0 



(87) 



Wc then have 



1 d 2 5r 1 d 2 Sr d_ 
~Z 2 ~d1?~ + z'dz 2 ' + dr~ 



Z d fr8r 
r dr \ Z 



1 r 3 

Z¥ \ dz 



2. Non-static spacetime (g t ^, ^ 0) 

When the spacetime rotates (gtcf, ^ 0), one finds that the Alfven wave ansatz (|81|l fails to satisfy the force-free 
equation to linear order in e. Physically, this is because the torsional and radial modes of the magnetofluid are 
coupled through conservation of angular momentum. One must therefore introduce simultaneous perturbations in Sr, 
S<f>, and 5t. In terms of these variables, the field perturbations become (to first order) 

SF™ -— (F<$8rY 5F W - ^ll F {0) ■ SF {1) - &-F<V ■ (89) 



9 4)st = 5^(^)st. (91) 



Here 

\9 J \9^J 

In this situation, the t-(f> components of the metric are not diagonal: g u = g^J det(t, 0), g^ = gtt/ det(t,(f>), and 
g t4> = -gt^l det(t, <t>), where det(i, <j>) = gttg<w - g 2 4 - 
The force- free equations are, again to linear order, 

J r{1) F { $ =0= J r(1) i^ 0) , (92) 

or equivalently 

J r(1) = 0; (93) 
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and 



Finally, the currents are given by 

4^V^/ t(1) 



d r (V=dg tt 9 rr SF t y 



(94) 



d z (V^9 tt g zz SF^) + d z (^g^g zz SF^); 



(95) 



^^gJ r{1) = QtiyTgg^SF^ + ^gg rr g u SF^) + d z (, 



-gg rr g z *5F^); 



(96) 



and 



47T A 



jj^ = d r {^—gg^g rr 5F^ + yf^g^SF^) - d t (^(g^) 2 5F^) 



Substituting for the field perturbations, Eq. I|93[) becomes 



1 dt 2 



= _ tt d_ (g^ dSr 
9 dz 2 9 dr \g^J dt ' 



(97) 



(98) 



One quickly obtains an integral of motion from this equation when the gradient in the z-direction vanishes: 



dt dr yg^, 



(99) 



The torsional perturbation 8<f> is slaved to the compressive perturbation dr through conservation of angular momentum. 
The manipulation of the force- free equation (|94ll is a bit more complicated; some details are given in Appendix B. 
The final result is 



, d 2 Sr 

lit? 



, d 2 Sr 
dz 2 



F, 



(o) 



+ 



d 



-gdr 



9 -^{5rF^ 



g<p<pg r r dr 



dt dr \g^ 



9_ ( gt±) F (o)_ 



dr\g++J r * 



(100) 



The wave operator on the left-hand side of Eqs. (|98|l and (|100|l is now mixed with an additional term involving 
the other coordinate perturbation. This source term can be treated as a perturbation when the spacetime is slowly 
rotating (g 2 , <C \gttg<t><p\)- It is also worth noting a simplification which takes place when the perturbations are uniform 
in z: then the right-hand side of Eq. I|100[) vanishes and, effectively, the fast mode equation involves only a single 
coordinate perturbation. That is, the time-evolution of the angular perturbation 5(f> is fixed by conservation of angular 
momentum, and does not react back on the equation of motion for Sr. More generally this is not the case, and the 
dynamics involves two independent fields. 



B. Scattering of a cylindrical Alfven wave in a spherical gravitational field 

A uniform magnetofluid is perturbed by a spherical gravitational field fSection lll D(l . We now consider the pertur- 
bation to a torsional Alfven wave propagating along this background magnetic field. We focus on a weak, spherical 
gravitational field, and allow for the possibility that the gravitating mass is extended in radius. Then the metric 
departs from the Schwarzschild metric; we write it as 

-gtt = 1 + A tt ; g rr = 1 + A rr . (101) 

(In this section, as in Section Til Dl r is the spherical radius.) In the weak field limit (|A tt |, |A rr | <C 1) the radial 
coordinate of the magnetofluid is perturbed by an amount Sr/r = 0(A). 
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The static perturbation 5r is a solution to the linearized version of Eq. I|27|) , 

d 2 (rSr) n (rSr) 1.. A . 1 d , , . . 



(102) 



Notice that the right-hand side of this equation vanishes in the Schwarzschild metric, A rr — —Au — 2GM/r. When 
the enclosed gravitating mass M is itself a function of radius, Eq. 11021) becomes 

d^_ 2 ( I 6r 1 = 2G dM 
dr A r z dr 

We expand the background magnetic field F(°\ the Alfven wave 5F^\ and its current SJ^ in powers of A, 



F (0) =F W a+f ( 



(0) 



-A 



and 



5F& = SF™ A + SF^ A ■ <5J« = s4 l J A + 5J[ 



-A' 



so that iv_A, an d SJY^* are first order in Sr. For example, expanding 1Z 

the background field, we have 



F, 



(o) 



0-A 



So r sin 0; 



F 



(0) 



0-A 



and 



F 



(0) 



l-A 



d(r5r) . 2 n 
= B -^j— ^ sin 2 0; 



F 



(o) 

l-A 



-Bo J" sin cos 9, 



2Bq rSr sin cos 8. 



(104) 
(105) 

2r(5r in Eq. for 
(106) 

(107) 



In the absence of gravity, the torsional wave propagates along the magnetic field. The associated charge and current 
densities 



5p = SJq(z — t) exp 

5 J = 5 J Q (z — t) exp 
are localized within a cylindrical radius rsin# ~ i?o- 



1/7- sin 
2 



i?0 

1 / r sin ( 
2 



l?o 



(108) 



We will focus on the case where the wave is strongly sheared, i.e., where its frequency is small compared with R 



Ro 



5J' Q {z-t) 



< 1. 



SJ (z-t) 

To lowest order in the gravitational potential, the components of the fluctuating current and field are 



5 7 {1) 

OJ 0-A 



= 5J\ 



5 J, 



(i) 

0-A 



cos 8 J] 



51 



(i) 

0-A 



= — - sin 95 J; 
r 



and 



SF, 



(i) 



0-A 



4tt5 Jo Rl 



Or 



sm f 



1 — exp 



1 / r sin ( 
2 



5F 



(i) 



0-A 



SF, 



(i) 



0-A 



Or 



SF, 



(i) 



0-A 



10 



SF 



(i) 



0-A 



Or 



(109) 



(110) 



(111) 



The static gravitational perturbation to the background magnetic field is known implicitly through the solution to 
Eq. (|102|) . We next solve for the perturbations SF^]} A and 5j[_ A . The current perturbation is determined from the 
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linear component of the force-free equation. Making use of 
are second-order in A, we have 



5 J, 



(i) 

0-A 



F 



(0) 



0-A 



fl IS 



= and neglecting the terms that 



OJ 0-A 



F, 



(0) 



l-A 



SJ 



(1) 

l-A 



F 



(o) 



0-A 



= 0. 



The force-free equation (|112Jl may be re-expressed by substituting Eqs. I)l()6[l . I|107(l and l)ll()[l . 



SJ 



(i) 

l-A 



SJ 



(i) 

l-A 



• cos 9 = —SJ(r sin 8,z,t) 



d{Sr/r) 



dr 



sin v cos f 



In the low frequency regime, the equation of current conservation simplifies to 



sm t 



5J 



(i) 

l-A 



= 0. 



(112) 



(113) 



(114) 



sin 6 89 

For the particular choice (|l(J8f) of incident current density, the solution to the combined equations (|113fl and l|114l) is 



(i) 



and 



SJ\ 



8J {1) 

0J 1-A 



Sr _ rSr , 9 

2 — cos 7T- cos sm" ( 

r B* 



1 d(rSr) Sr , o 



<5J 



(115) 



(116) 



Notice that the current SJ^]} A receives a contribution from the metric perturbation as well as from SF^ A , e.g. 



47T 



8I {1) 

0J 1-A 



1 d 



{smeW 1) M-aI } + -5^-7 |j {sine [*f£? a 1 } . (117) 



Here, the metric perturbation is g rr — 1 = —2GM/r to lowest order. Thus the current components [<5</[l^] r ' e are 
sourced by the fluctuating field 



SF 



(i) 



l-A 



= — \A-K(r5r) sin 0] <5 J(r sin 6, z,t) + 



2GM 



SF 



(i) 



0-A 



(118) 



The gravitational perturbation of the charge density, and the electric field which it sources, can be obtained by 
combining the Bianchi identity 



do 



6F£ 



(i) 



dr 



SF, 



with the MHD condition 



5F {1) 

Or 0-A 



10 



SF 



(i) 

l-A 



l-A 



t0' 



SF 



0-A 



SF 



(i) 

l-A 



10 



F 



F 



di 



(0) 
0-A 

(0) 
0-A 



SF^ A 



Or 



<j>r 



F, 



AO) 

l-A 

(0) ' 
l-A 



0. 



(119) 



At leading order in the gravitational potential, this last equation becomes 

(i) 1 \ SF w> l —A 1 ** 1 -) 



SF 



l-A 



= r cot ( 



10 



l-A 



— cos L 
tr d In r 



SF 



(i) 

0-A 



Or 



Substituting this expression into Eq. (|119fl . one can solve for the perturbation to the electric field, 



SF, 



l-A 



-47r<5J 



F() 



d(Sr/r) 
dr 



(l-e- R2 ^ R o)+ S in 2 eSre- R2 ^ R o 



and 



SF 



(i) 

l-A 



10 



= -[4n(r5r) sin9cos9]5J e- R2/2I %. 



(120) 



(121) 



(122) 



(123) 
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Here R = r sin 9. The charge density perturbation is then 

r 2 dr 



47T 



<5/ (1) 

OJl_A 



[r- 



5R 



(i) 



l-A 



+ 



1 



1 6*96* 



sin 9 



l-A 



{- 

{sin9(g u + l) 



} 



SF 



(i) 



0-A 



10 



r 2 sin 6* 96> 

Here <7** + 1 = —2GM/r to lowest order. 

Our principal goal here is to work out the second-order current [8J?^ A^. This may be obtained from 



(124) 



SJ, 



(i) 

0-A 



5F, 



(i) 



l-A 



SJ 



(i) 

l-A 



SF 



(i) 



0-A 



SJ 



(2) 
l-A 



F 



(0) 



0-A 



(125) 



by substituting Eqs. (110), (111), (116), (118), (122) and (124). The resulting expression is fairly complicated, but it 
simplifies dramatically when the gravitating mass is concentrated at a radius much smaller than the transverse scale 
i?o of the incident Alfven wave. In this regime, one has Sr/r ~ G(dM/dr) <C GM/r, and the dominant field and 
current perturbations are 



2GM r 



Or 



SF, 



(i) 



0-A 



/ 2GM \ AnRl 
\ r J sin 9 



1 - e - R2/2R o 



SJo(rcos9 — t), 



and 



SJ 



(i) 

l-A 



2GM 



co S 2 9e- R2 ' 2R --^-{\ 



_ e -R 2 l1Ri 



SJo(r cos 9 — t). 



Combining the terms on the left-hand side of Eq. I|125|) , and substituting the background field H19|) , we find 



SJ {2) 



iTr(SJ Q ) 2 (GM\ (R 



Bo sin 2 9 V r 



-R 2 /2R 2 



R 2 /R 2 



,R 2 /2Rl 



This expression simplifies further at small cylindrical radius (R <C Rq) 

\6J \ 2 (GM\ 



SJ 



(2) 
l-A 



6tt- 



Bo \ r 



sin 2 6 e 2 ^"*), 



(126) 



(127) 



(128) 



(129) 



when the incident Alfven wave is a pure harmonic. The second-order current is the source for an outgoing compressive 
wave 



d 2 {r5n) 
dt 2 



1 



n2 sin 9 d8 



d ( d(r5rj)\ d 2 (r8rj) 
sin 9 



09 



dr 2 



4nr 



2 r 



SJ 



(2) 
l-A 



(130) 



of amplitude Sri . As in the case of a torsional wave propagating in a cylindrical spacetime, the compressive mode is 
second harmonic and second-order in the wave amplitude. 

The Greens function solution to Eq. i|130|) is dominated by the fluid at radius r ~ R . As a result, the shape of 
the asymptotic fast wave is sensitive to the distribution of current within the incoming Alfven mode. We will not 
examine the pattern of the outgoing compressive wave in full detail here. 



V. COLLISIONS OF AXISYMMETRIC WAVES IN A CYLINDRICAL SPACETIME 

In this section, we study the collision of two torsional Alfven waves in a uniform magnetofluid in a static cylindrical 
spacetime (g t $ — 0). The effects of spacetime curvature are subtle enough that we specialize to axisymmetric modes, 
supported only by a perturbation S<p. The result easily generalizes to the case of flat space, and in the next section 
we examine more general types of wave interactions in a Minkowski background. 

The two Alfven waves, propagating oppositely along the magnetic field, are labeled + and — , and the component of 
the current generated by their interaction is labelled /. The colliding waves are each assumed to be supported within 
a cylindrical shell, between radii i? m i n and i? max . We start from the force- free equation 

( J\ + Jt) (SF+ + SF t - ) + (J* + Jt ) (SF+ + SF-) + jfF$ =0. (131) 
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The background magnetic field is defined by J^(°) = [Eq. Q32[)]. The fluctuating fields are 



>± F (o). 



92 0r ' 

r flir \~"z f J r zr , 



0f 



where t> 2 is given by Eq. (JSUJ. The associated currents are 

1 d 



4:TTy/—g dr 



-gg tt g rr 5F±) = ± 



1 



47T\/— g dr 



-SFl 



and 



1 



d 



Airyt—g dr 

Combining these expressions, the force-free equation becomes 

4?r jj, _ 9rr_ d 2 Sr grr d 2 5r ^ d_ 
B Q V 9 1 ~~ gu dt 2 g zz dz 2 



-gg zz g rr 6F^ 



dr 



g_d_ ( Sr 9rrg<t>4> 



d_ 

9<p<p) dr 



fjz 



9zz 



t d(t) + dcj>- 
dz dz 



dz dz 

2 



dz 



d<\>- 
dz 



d r (In v z ) . 



(132) 



(133) 



(134) 



(135) 



In contrast with the case of a single (unidirectional) Alfven wave, a source term appears even when the wave speed 
has a vanishing gradient (d r v z = 0). For that reason we must, in general, retain the curvature corrections in the wave 
operator for Sr. 



A. Collision in a nearly flat, cylindrically symmetric spacetime 

We now specialize to a nearly flat cylindrical spacetime with Newtonian potential A and metric l|34|) . and focus 
on the source terms in Eq. I|135|l that are proportional to {d z <f>+)(d z (f>-). The torsional wave speed is v z = 1 + 2A. 
Comparing with flat space, the new effects wc arc interested in arc already present in the case of a very long wavelength 
twist in the magnetic field, 



dz 



«1, 



so we neglect the derivatives in t and z. The perturbation to B z is 



SB Z 



SF^ r _ l d{r5r) 



? (o) 

(fir 



and the equation of motion i|135D simplifies to 



d_ 

dr 



SB Z 



2_d_ 

r 2 dr 



r dr 



dz dz 



- 2— dr, 
dr 



4r z 



(136) 



(137) 



(138) 



dr dz dz 

Our first goal is to calculate the amplitude 6R max (t) = Sr(R max ,t) of the hydrodynamic fluctuation at the outer 
boundary of the annulus supporting the Alfven waves. The amplitude of the fast mode at large radius can then easily 
be obtained, to leading order in A, from the equation of motion (|61|l in Minkowski space. Taking the two waveforms 
to be [3 



-.(t,r, z) 



t> ± sin 



k(zTt) f±±(r), 



(139) 
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one sees that the interaction term 

d(f) + d(j>- 
dz dz 



1 



-k S<p 0+ 54>o- 



cos(2kt) + cos(2kz) f±+(r)f±-(r) 



(140) 



is the sum of a uniform oscillation of the cylinder, and a static deformation that excites no wave motion outside the 
cylinder. Thus the calculation of the fast mode amplitude reduces to the problem of calculating the amplitude 5R maK 
of the oscillation at the outer boundary of the cylinder. 

We start by integrating Eq. Q137|l outward from SR = at r = to obtain 



mm -R 



^-r 2 dr 
dr 



SB Z 



(141) 



to first order in A. The expansion of the magnetofluid is nearly incompressible in the near zone outside the cylinder 
(kr <C 1), so the relevant inner and outer boundary conditions are 



0. 



and 



dz 



= 0. 



The field perturbation is nearly uniform inside radius i? m in; we will only need the leading term 



6B Z 
"Bo" 



-Rmin 



l_d_ 

r 2 dr 



dz dz 



dr = 4$ + _(i? max ), 



where 



$ + _(r)= f : 



-dr. 



dz dz 

The radial perturbation within the annulus i? m i n < r < i? ma x can be found iteratively from 

~6B Z _dA 
"Bo" 



(142) 



(143) 



(144) 



(145) 



rSr 



-^min^-^mni 



Rn 



2-5r(r>) 



dr' 



1 2 SB Z 
2 r ~Bo" 



1 -r^( S ^)dr' 



2 dr' V Bi 



r' — 5r(r')dr' , 
dr 



(146) 



after substituting Eq. I|138|l in the third term on the right-hand side. Evaluating this expression at r = i? ma x we find 



RmaxSR 



max UIt max 



,d\ 
dr 



2* + _(r) 



dr 



dr + 4$ + _(i? max ) 



2 d\ 

r —dr. 

dr 



(147) 



The first remark to make about this expression is that <5i? max vanishes in flat space (dX/dr — constant). The same 
turns out to be true in the special case of a cylindrical line mass (dX/dr = K/r), but not for a more general cylindrical 
mass distribution that is extended in radius. 

We now obtain the fast mode amplitude in the wave zone fer> 1. The cylindrical wave equation for SB Z has the 
outgoing wave solution 



^ = Ae- tkt H {1) (kr) ~ A\ J_ e «=0--t)-W4 f kr » i) 
Bq V Trkr 

r(i) 



(148) 



in Minkowski space. The Hankel function scales as Hq (kr) ~ (2i/tr) ln(fcr) at small radius. To obtain the normal- 
ization factor, we note that in the near zone (kr <C 1) where the magnetofluid is nearly incompressible, the rate of 
transport of magnetic flux is nearly constant in radius, 



dt dt 



constant (kr <1). 



(149) 
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This implies 



dSr 
~dt 



d{5r) 



dt 



Substituting this expression into E — —dt(5r)Bo(f x I), and thence into 



and integrating over radius, gives 

SB Z SB z (R n 



and 



The general solution is 



Bo 



d 2 (Sr) 



dt 2 



In 



Ae~ tkt = y (^ 2 i? max ) Sr(R max ,t). 



SB z (r,t) 

Bn 



t-r o2 



d 2 (Sr) dt' 



^max i 



dt 12 ^{t - t') 2 - r 2 



(150) 



(151) 



(152) 



(153) 



(154) 



B. Collision in a static black string spacetime 



As another illustrative example, we consider the collision between two torsional waves in the spacetime of a static 
black string [Eq. with J = 0]. The background magnetic field F^J = B Q r/(r 2 /£ 2 - M) is aligned with the 

axis of the string. We assume, as before, that the colliding waves are supported only within a cylindrical annulus 

-Rmin < 7" < R max , and take i? m in, Rmax ~ £• 

A radial disturbance of frequency id <§; t~ x accumulates a phase 



1/2 

| dr = u 

9tt 



fJ 



dr 



r 2 /£ 2 



M 



(155) 



This integral converges at large radius but diverges near the horizon rn = M x /H. When M < (uj£) 2 , there is a well 
defined wave zone which is localized at a small radius, r ~ u>£ 2 . The effect of a small finite angular momentum J is 
to impose a reflecting barrier for the magnetosonic wave at a radius r ~ (J 2 /4w) 1 / 3 , and therefore to form a resonant 
cavity. We will not consider the case of finite angular momentum here. 

The medium outside radius i? m ax will, at the same time, respond to a low-frequency disturbance with a uniform 
compression or rarefaction, SF^r ~ F^. . Thus Sr ~ r at a large radius (r ^> £). The amplitude of the radial 
perturbation SR max can be obtained by noting that the magnetic flux 



$ = / 2irdrF r<f) = / dr 



2nrB t 



r 2 /£ 2 -M 

diverges at large radius. At the same time, the background field energy 



dE B 

dz 



rrQ(f>(f> 



F 2 

7 " „rr <j><t>_L± 

J y y 8tt 



rdr 



4 {r 2 /£ 2 -Mf/ 2 



(156) 



(157) 



of the outer boundary of the annulus is accompanied by a small field 



is finite. As a result, an oscillation 5R„ 
perturbation, 5B z (R max ) ~ 0. 

The nature of the wave solutions at small radius is most easily demonstrated by taking the limit M — * 0, Z(r) 
r 2 /£ 2 . Then the wave equation (|75|) . which applies to the case of a uniform radial oscillation Sr(R max ,z,t) oc e~ u 
simplifies considerably. It transforms under the change of variable r* = — £ 2 jr to 



d 2 (6B Z 

~di? Kb 1 



_d_ 

9r» 



1 d 



8B Z 
~B Z 



= 0. 



(158) 
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A further change of variable SB Z /B Z = h{r*){u>r*)e %u}t yields the Bessel equation 

nd 2 h dh 



r: M + r *^ + ( " 2r *- 1)/l = ' 



which has the wave solution 

SB 2 
IF 



7r(wr* 



-Kj(r,+t)-37r/4 



(159) 



(160) 



propagating to small r (large negative r*). It will be noted that the amplitude of the wave diverges in the wave zone, 
because the radial phase velocity 



k 



9tt 



(161) 



asymptotes to zero. 

To relate the coefficient A to the amplitude of the oscillation at radius i? m i n , we note that the rate of transport of 
magnetic flux 



2nF, 



d{Sr) 2irB £ 2 d(6r) 



dt r dt 

is constant in the near zone (|wr*| <C 1). Hence d r (Sr) ~ 8r/r. From Eq. (73) we have, 

d 2 5r* d f5B z 



dt 2 dr* V B z 

Combining this with the low frequency expansion (u;r*)/fi (a;r*) ~ 2i/iz — (i/n)(ujr*) 2 ln(o;r*/2) gives 

^ -<7r/2 8R* min (t) 



Ae 



(162) 



(163) 



(164) 



The last step in obtaining the amplitude of the fast mode that emerges from the wave collision, is to express SR mul 
in terms of the amplitudes of the colliding Alfven modes. The terms involving (d z 4>+)(d z 4>-) m Eq. IJ135JI can be 
combined to give 



d fSB z 



d_ 

dr 



d ( Sr 



Integrating with respect to radius, this becomes 



SR„ 



5R„ 



Ru 



Rn 



= 2 



dr 



drr 



dr \ r 
dr 

1 - In 



2 d 



r dr V dz dz 



(165) 



l -Rmin ) 


SB 2 ' 


-In 




~B* _ 


in ) 


1 d^+ 


d<t>- 




dz 


dz ' 



Ru 



d (SB 2 
dr VB 2 



(166) 



since d z <j)±{R m \ n ) = <9 z 0±(i? max ) = and SB z (R max ) = 0. The field perturbation at the inner boundary is obtained 
by substituting i? m i n — > i? m ax in the logarithms in the first line of Eq. I|166(l . After subtracting the two equations, 
one obtains 



SB 2 
"B 2 



-Rmin 



d<f> + d<j)_ 
dz dz 



dr. 



(1666) 



VI. COLLISIONS OF NON-AXISYMMETRIC WAVES IN MINKOWSKI SPACE 



We will now broaden our discussion of mode collisions to non-axisymmetric modes, and consider both fast and 
Alfven modes. Three-mode interactions in the force-free limit were previously considered in |3j, using the magnetic 
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Lagrangian formalism developed there. We return to this subject and clarify the properties of i) the three-mode 
coupling between two Alfven waves; ii) the interaction between two fast waves; and iii) the interaction between a fast 
wave and an Alfven wave. 

Basic constraints on the mode interactions arise from conservation of energy and momentum. Because the relativistic 
magnetofluid is inherently compressible, there is guaranteed to be a three-mode coupling between two Alfven waves 
and a fast wave |3j. The Alfven waves (Al, A2) satisfy the dispersion relation lua = ±k A ; and the fast mode uj f — \k F \. 
Conservation of energy 

&A1 + ^A2 = u F (167) 

and longitudinal momentum 

k Al + k A2 = ^Al ~ UJA2 



k z F =oj F cos6 F (168) 



(cos8 F — k F ■ B /k F B ) can be combined to give 



U)A1 — UA2 ncn\ 

cosOp — . (169) 

UJAl + UJA2 

Thus in a frame in which the two colliding Alfven modes have equal frequencies, the fast mode is emitted perpendicular 
to the background field and with a frequency uj f — 2lua ~ just as we found in Section IV Al 

There is also a three-mode interaction between the two colliding Alfven waves and a third Alfven wave. From a 
kinematic view p o int, this interaction is non- vanishing only if at least one of the colliding modes has a zero-frequency 
component ([14(, [111), as is easily checked by setting cos9f = 1 in Eqs. I|167|) and i|168|) . An alternative description is 
that the magnetic field lines experience a net displacement (or braiding) across at least one of the colliding wavepackets 
|16| , Q • (A further requirement for a non- vanishing three-mode interaction, as we detail shortly, is that the colliding 
modes are non-axisymmetric.) 

Next consider the interaction between two fast waves. It is here that a significant difference arises between these 
two viewpoints. If the background magnetic field suffers a net displacement across the fast wavefront, then a zero- 
frequency component is present in the fourier decomposition of the field. From a kinematic viewpoint, this zero- 
frequency component would facilitate a three-wave interaction between the two colliding fast modes and a third fast 
mode. Two colliding fast waves of finite frequency can generate a third fast wave only if the waves are colinear: only 
in that case is it possible to satisfy conservation of energy 

Ul F l + LJ F2 = UJ F3 (170) 

and momentum 

kfi+k F2 = k F3 . (171) 

(Two photons propagating obliquely in vacuum cannot merge to form a single photon, in part because these kinematic 
conditions cannot be satisfied.) From the same kinematic viewpoint a three- wave interaction between two colliding 
fast waves and an Alfven mode is possible only in the degenerate case where both colliding waves propagate along 
the background magnetic field. This restriction would appear to be somewhat relaxed if one of the waves has a 
zero-frequency component: then the direction of propagation of that wave can be arbitrary. 

As our second task in this section, we calculate the form of this three-mode interaction for two colliding planar 
fast waves. We find that, in fact, no fast wave is emitted during the collision, but an Alfven wave is. The Alfven 
wave is emitted for any direction of propagation of the two colliding waves. Moreover, when there is a field line 
displacement in one fast wave, and the other is a pure sinusoid, then the outgoing Alfven wave is also a pure sinusoid 
but its frequency is not obtainable from a three-wave resonance condition. It is only when the two colliding modes 
are directed along the background magnetic field that the resonance condition can be satisfied with a zero-frequency 
component. 

The fast mode also undergoes a four-mode interaction analogous to photon scattering, with the conservation relations 



u; F1 + uj F2 = u} F 3 + u F4 (172) 

and 

k F i + kp 2 = k F 3 + k F 4- (173) 

However, the fast mode also has a non-linear interaction with an Alfven mode. We show that if the Alfven wave 
spectrum is that of a critically balanced cascade (|17|), then this second interaction dominates the self-interaction of 
the fast mode (as was claimed without detailed justification in Q). 
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A. Collision between two Alfven waves 



A collision between two axisymmetric torsional modes generates a compressive motion transverse to the axis of 
the background magnetic field. The amplitude of this compressive motion is, however, strongly suppressed when the 
colliding wavepackets are highly elongated {k z R max <C 1). In the case of two monochromatic torsional waves with 
frequencies oj = ^k Z: this radial disturbance vanishes at the outer boundary of the zone supporting the colliding 
Alfven waves. More generally, <5i? max is suppressed by a factor ~ (k z R max ) 2 , compared with the result of simple 
dimensional analysis, 



5Rn 



R„ 



{k z R n 



) 2 xi?i 



2 d(j) + d4>- 
max dz dz 



(174) 



In other words, the three-mode interaction between two Alfven waves and the fast mode is suppressed for elongated 
wavepackets (such as are created during a weak turbulent cascade 0D- 

There is, nonetheless, a much stronger three-mode interaction involving a third Alfven wave. This interaction is 
not apparent if one restricts the two colliding waves to be axially symmetric, and so we will broaden the analysis here 
to include Alfven modes supported by coordinate perturbations 8cj> and Sr. The relevant component of the force-free 
equation is 



7*(i) F (1) + 7 Z W F {1) + J r ^ F {0) - 



(175) 



In Minkowski space, 5F^ — (Bor)dtSr and SF^' — (Bor)d z Sr. The second-order current vanishes for the two 
separate modes + and — , but where they overlap the first two terms in Eq. 1)175(1 do not cancel: 



J*MF$ + J Z[1) F^ = 5> [d 4 (d z Sr + d z Sr^) - d z Sr+ d r (r 2 d z 8cj>^ - 8 z Sr_ d r {r 2 d z 8cjy + )] . 

When at least one of the colliding Alfven wavepackets is not axisymmetric (with both perturbations 
excited), one also finds an explicit second order contribution to the fields, 



(176) 
and 5r 



—d r (rSr + ) d t S(j)- +-d r (r5r^) d t 5<p + — d t Sr + d r 
r r 



dtSr^d,, 



r"<P+ 



Bnr; 



(177) 



(2) (2) 

and similarly for SF^J and SF 7 - Z ' . The second-order current therefore has an explicit contribution from the colliding 
modes. As before, an additional interaction component Srj, 6<f>i to the coordinate fields is required to solve the 
force-free equation: 



4nJ r ^ = -dtSF™ +d s 5F<? + ^5F { r 



(2) 



-2B r 



d z Sr + d r d z S(j)- + d z 6r-drd z 5(j)+ + (d z 



'->+U z i 



B r <j -dfS^i + d 2 z t 



^d r (rSn) + 1 



(2) 
rcf> 



(178) 



Thus Eqs. 1(175(1 and (1 176(1 combine to give the equation of motion 



d 2 S<f> I + d 2 



2 
r 



-d r (rSn) + d+Sfa + -^sf^ 

r Bqt 



(2) 



2d z Sr + d r (rd z 5(f)-) + 2d z 5r- d r (rd z 8<f) + ) + rd^, (d z 5(f> + d z 



-) ^ (d z 5r + d z 5r-) 

r 



(179) 



In this equation, the transverse components of the laplacian on the left-hand side suppress the interaction unless 
they vanish - that is, unless the the new mode is an Alfven wave (Eq. I47|) . Even without making this restriction, the 
dominant three-mode interaction can be obtained by transforming to the variable 



T± = d> d(j> 5<f>±((j>,r)] 



(180) 
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Here, the integral is performed at constant cylindrical radius. Applying the operator § dip to Eq. \ 179(1 kills off the 
last term on each side, and we are left with 

-<9 t 2 r/ + d*Tj = -j> d(j> [d z Sr+ d r {rd z 5<j)-) + d z 8r_ d r (rd z 5cj> + )} . (181) 

Transforming to light-cone variables z± = z ± t, d± = \{d z ± dt), the three-mode correction to each Alfven mode can 
be calculated from 



0_r+| out -a_r + | in = J ' dz+d+d-T! 



= i jdxj> [d z Sr+ d r (r<9 2 A0_) + d z Ar_ d r (rd z S(p+)} . (182) 

The correction to T_ is obtained by interchanging + and — in this equation. 

As advertised, the resonant three-mode interaction depends on the presence of a net twist (radial displacement) of 
the magnctofluid across the wavepacket, 

A0± = / dz T d z M ± (z T ); Ar± = / dz T 8 z 5r ± {z T ). (183) 



Because the colliding modes are Alfven modes, Ar± is related to A<p± through the condition of incompressibility, 

^d r (rAr±) + <9 A</> ± = 0. (184) 

The three- mode interaction also depends essentially on the non-axisymmetry of the colliding wavepackets (i.e., on the 
simultaneous presence of coordinate perturbations dr and 54>). Finally we note that applying the operator J dip to 
the constraint equation (|48|) for the fast mode gives 

^ (V 2 J dcpSc^j = 0, (185) 

hence r 2 r = constant. As long as T does not have a singularity at r = 0, we see that r = for the fast mode. 



B. Collision between two fast waves 



In order to tackle the collision between two fast waves, we consider the special case of two planar waves. (There are 
inherent complications associated with the choice of cylindrical geometry, due to the the non-locality of waveforms 
propagating in two spatial dimensions - in contrast with one or three dimensions, the Greens function is not a delta 
function.) 

Choosing a background magnetic field Bo = Bqz, the perturbed Faraday tensor is 

F^u = (d^d^y - ^y°cU°) B . (186) 

Each wave, by itself, is equivalent to a plane electromagnetic wave superposed on a uniform background magnetic 
field, and the associated coordinate perturbation satisfies the usual wave equation, 

x° = x + Sx° li2 {x^); d v d v [Sx° 12 (x^)] = 0. (187) 

We can boost along the background magnetic field into a frame in which one of the colliding fast waves propagates 
in a direction perpendicular to Bo, say 

= 5xl(x-t)x. (188) 

A specific example of such a waveform is the harmonic perturbation ki = k\x, 5x\ = SXie lkl ^ x ^ tS> . Because neither 
fast mode involves a propagating current in the force- free limit, the interaction vanishes if the electric vectors of the 
two modes are parallel and their direct superposition involves no violation of the constraint <JE • Bo = 0. Thus the 
waveform of the second mode is chosen to have a non-vanishing derivative in the y-direction. (In general all three 
components of k2 are non- vanishing.) The corresponding coordinate fluctuation is 



5x°=<52/°(n 2 -x-t) 



y h x 

n 2y . 



(189) 
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where n 2 is an arbitrary unit vector. The relative normalization of the fluctuations in a; and y° is fixed by the 
vanishing of the current, d y Sx 2 — d x Sy 2 = 0. We allow for the possibility that there is a non-vanishing displacement 
across (at least) one of the colliding waveforms, e.g. 

Ax? = / di(8xl)'{t) ± 0; Ay° 2 = / dt{5y 2 )'{l) * 0. (190) 



The current obtained by superposing Sx° = <5x? + <5x 2 is 

4ttJ£ f = d v (d»x°d v y° - d' 1 y d u x ) B = {d^Sx^Syl - d u d"5yld v 8x\) B . (191) 

Because the linear currents J^ l > vanish, the force-free condition is J x ^ = J y ^ = to second order in the wave 
amplitude. We must then introduce an interaction field x^® and associated current 

Air J? = [dyd'x ! - d x d^y°j + 6% (d v d u Vl ) - 5$ (d v d v Xl )] B . (192) 

The conditions 

Jf + J FF =Q = Jf + J y FF (193) 

correspond to the equations of motion 

d v d v yj + d x (d y Sxj-d x 5yj) = 

(1 - n 2x ) [{8x\)"{x - t) ■ (<5y 2 °)'(n 2 ■ x - t) - n 2x {8x\)'{x - t) ■ (<Sy 2 °)"(n 2 -x - t)] 

(194) 

and 

d v d v x°j - d y (d y 5x j - d x 8yj) = n 2y (l - n 2x ) ■ ~ t) ■ (5y° 2 )"(n 2 ■ x - t). (195) 

It is useful to project these equations onto the outgoing fast (d x x°j + d y y T ) and Alfven (d x y°j — d x y T ) modes, 

d v d v (d x x°j + d y yj) = 2(1 - n 2x )n iy {8xl)" \x - t) ■ (Sy° 2 )"(n 2 ■ x - t) (196) 

and 

{-d 2 t + d 2 z ) (d x y°j - d yXl ) = (1 - n 2x ) f {8x\)'"(x - t) ■ (8y°)'(n 2 ■ x - t) - 

Ay) 



("L + n 2 2y ){8xl)'{x - t) ■ (<$y°)'"(n 2 ■ x- i) 



(197) 

To find the outgoing perturbation to the fast wave 1, we start with the Greens function solution to Eq. (|196fl . 

(d x x°j + d yyi ) (x, t) = 2n 2y (l - n 2x ) J d 3 x' --L-(fo°)" \x' - t + |x - x'|) • (<5y 2 °)"(n 2 • x' - t + |x - x'|), (198) 

(using the retarded time t — |x — x'|) and take the limit x, t — ► 00 at y = z = 0. Each wavepacket is localized 
within a distance ~ L along its direction of propagation. Transforming x' — > x' — x and then setting {x',y',z'} — 
r{sin 8 cos 4>, sin 8 sin (j>, cos 8} , one has 



(d^j + dyyf) (x,0,0,t) = 2n 2y {l - n 2x ) j rdr sin 0d6dcf> { (8xl)"[r(l + sin cos </>) + x - t] 

(Sy%)"[r(l + n 2 -n)+n 2x x-t\y 



(199) 



The integrand is non- vanishing only in a small region Ar ~ L near r = (t — n 2x x)/(l + n 2 -f2) and A6A<j> ~ L/t near 
8 = 7r/2, 4> = ir. When integrating over r, the argument of iSx®)" varies only by a small amount ~ L 2 /t. The integral 
J d£(Sy 2 )" (£) itself vanishes, and so the outgoing fast mode has an amplitude ~ (L/t)(8x ( ()'(Sy 2 y — > as t — > 00. We 
conclude that the two fast modes do not couple to a third fast mode, even if one of the colliding waveforms has a 
non- vanishing displacement l|19U|) . 
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Now let us examine the coupling to an outgoing Alfven mode. Transforming to the light-cone coordinates z± = z±t, 
and taking x = y = 0, Eq. (|197|l becomes 



+ & 



.(fl x y? - a^?)(0,0,2+,z_) 



0\/ 



(Sy Q 2 Y 



\( n 2z + + ^(«2 2 - 1)Z+ 

1. , 1. 

-(n 2z + 1)2- + -(n 2z - l)z+ 



(200) 



For example, the outgoing Alfven wave propagating to large positive z is 
cL(a,y? - d y xfj (0,0, z-t) = 

(hi) 1 



-(l-n 2x ) I dz + <!(&?)'" 



2(*-"*f) 



1. si, 

2 ("2a + l)^- + -(Tl2 2 - 1)2+ 



{nl x + nl v ){8xl)' 



2 (z -" Z+) 



(^2 



1, 1. 

-{n 2z + l)z_ + -(n 2 * - 1)2+ 



(201) 



This integral has a closed form solution when one of the wavepackets (say 2) is much shorter than the other and has 
a zero-frequency component. Then we can treat (Sy 2 )' as a delta function, (Sy 2 )'(£) = AyJM^). After restoring the 
dependence on the transverse coordinates x and y, and integrating once over z_, our final expression for the outgoing 
wave is 



(d^ - d v xf) (x, y,z-t)= { -K^-Ay° ■ (6x1)" 



2n 2 



1 - n 2 x 



n 2z 



1 - n 2z 



n 2y 

1 - n 2z 



n 2z 



n 2z 



-(z-t) 



The outgoing Alfven wave is sheared, with wavevector proportional to 

(w, k) oc (n 2z ,n 2x + n 2z - 1, n 2y , n 2z ) . 



(202) 



(203) 



Notice also that the amplitude of the outgoing wave involves a higher-order derivative of the ingoing wave, and so the 
higher-wavenumber fourier components are enhanced by the collision. 
The frequency of the outgoing Alfven wave is 



n 2z - nig 



n 2z 



1 



-U> F i 



(204) 



when the fast wave 1 is a pure fourier mode with frequency u)pi, and its wavevector has some component along the 
background magnetic field. The three-wave resonance condition 



U>A — Wfi + LUp 2 = LJfi] 



k F i 



l F2 



k 



(205) 



can be satisfied with u>f 2 = k F2 = if the propagation of fast wave 1 is aligned with the background field (ni z = 1). 
But, more generally, Eq. (|204Jl is not consistent with such a three-wave resonance: the integral over z + in Eq. 1201|) 
is non- vanishing only if 



(ni z - 1)wfi + (ni z ~ 1)wf2 = 0. 



(206) 



A solution to this equation, with at least one non- vanishing frequency, requires that either n\ z = 1 or n 2z = 1. 
This illustrates why a net field line displacement across an MHD wave has only a limited description in terms of a 
zero-frequency component of the wave. 



C. Collision between a fast wave and an Alfven wave 



The collision between a fast wave and an Alfven wave can generate both an outgoing Alfven mode and an outgoing 
fast mode through three- wave interactions. It is useful to compare the strength of this interaction with the interaction 
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between two fast modes discussed in the previous section. We focus on the particular case where the Alfven wave 
spectrum is highly anisotropic, and the coupling parameter k y x° — k x y° ~ (k±/ui)(5BA >u /B) is independent of scale 
and close to unity. A key difference between these two types of wave collisions is that the sheared Alfven waves are 
current-carrying, ~ k±SBA,ui and there is a second-order Lorentz force involving the interaction between this 
current and the fluctuating field of the fast wave. The magnitude of this Lorentz force is 

r(i) 



,m ~ k^8BA,uj5Bp L 



(207) 



Here 5Ba,u and 8Bp^ are the field perturbations of the Alfven and fast waves at frequency ui. We can compare this 
with the second order current generated in the fast-wave collision, 



B 



The ratio is 



7(2) r> 
■Jpp&O 



kj_SB AiUJ 

UjBn 



/SB 



F,u 



V Bo 



SB 



Bn 



> 1. 



(208) 



(209) 



Thus, if the Alfven waves form an anisotropic cascade, then the fast wave spectrum evolves predominantly through 
collisions with Alfven waves and not through self-collisions between fast waves. 
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APPENDIX A 

In this appendix, we show how the force-free equation may be derived from the action 

S' = I d+xoL'^ [ d A x-^T~gg w g va F^F"^ I d^x^^g^g^F^F"" , 



using the magnetic lagrangian fluid variables (|1(J|> . Varying Eq. IjAljl . 



5S' = J d i xo5L l = J d A xo 



dL' 



5x" 



dL' 



dx^ 



_dx» d(dx^/x Q a ) \dx a a 

there are terms which arise from the explicit x^-dependence of the metric, 



dL' 
'dx' 7 



(-go) dx 13 dx~* dx 5 dx e ~ A ~ vC 
4J 4 ftcg dx* 3x1 dxi ° ° 



d 


gpsg^e 


dx** 


. V~9 . 



as well as from the dependence of F^ v on dx^ 1 /di 



1 



dL 1 



d(dx»/xo a ) 



9l38g~<e 



-g 



d 



d(dx^/x a ) 



(-g ) dx 13 dap dx s dx e ~ kX ~ vC 
4J 4 dx$ dx$ da% dx c ° ° 



After extremizing, the usual Eulcr-Lagrange equations are obtained, 



d 
dxQ 



dL' 



d(dx»/dx%) 



dL' 

dx» 



Substituting eqs. (|A3|) and (|A4|) then gives 



d 
drt 



dx' 



~ va \ _ dJ 4 
o J d(dxf>/dx%) 4 



9 



-gpsgje 



pPnp< 



5e 



(-g)Ji 



4 dx^ 



= 0. 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 
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Now the background field satisfies 



9 ( r^dx s ~ m 
dx% \y 9o dxf° 



(A7) 



Making use of this relation, and the identities 



9 



dJ 4 



d{dx»/dx%) 



= 0; J 4 



9J 4 



9 



9 



dx^ d(dx^ /dxy ) Oxq 



(A8) 



in Eq. (| A6|) gives 



F a 



dx e 



2 dx* 



= 0. 



(A9) 



To demonstrate the equivalence of this expression and the usual force free equation, one need only substitute 
F a p = —\\f—ge a p pa F ( " y in 



f 9 



(A10) 



After making use of the contraction 



(-<?) 



(All) 



Eq. (|T3|) is obtained. 



APPENDIX B 



In this appendix we list, for reference, the axisymmetric perturbations of a background magnetofluid satisfying 
J p = 0, using the electric Lagrangian fluid variables (©■ 
The field strengths are given by 



r / dr o d0 o 90 o 9r \ (0) f 9r dtp dt Q dr \ (0) 
" " ~ dr d(f> dr d(j>) r<t > I dr d<j> dr dcfr J rt ' 



(Bl) 



where the background fields F^ and are 



F, 



(o) _ B g ^^ffo , 



'-gaRo 



(B2) 



(B3) 



In a cylindrically symmetric spacetime with z = z, dcj)o/d4> = 1, and dr^/dcj) = dt /d<fi — 0, the perturbations 
6F„ V — Fn V — Ffi) have the form 



(B4) 



_ /9r 9<£ 9</> dr \ (0) / 9r 9i 9< 9r \ (0) p (o), pV 



(B5) 



_ ' 9r 0<fo 90o 9r \ (0) , , / 9r 9i 9t 9r \ (0) , . . 
SFrz ~ I 977^7 ~ 797 97 J M + I 97 97 _ 7977*77 I Frt M ' 



(B6) 
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and 



fdr d(j>o d<f>o dr \ (o) , N . f dr dt 
Substituting F^(ro) = — {g^ / 9q)F^ {fo), these equations reduce to 

IP _ 5r P (0)/ \ p(0)/ n 



dt dr \ (o). 



(B7) 
(B8) 

(B9) 
(BIO) 
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APPENDIX C 

This appendix is devoted to a derivation of the force-free equation (|100|) in a non-static, axially symmetric spacetime. 
The force-free equation (|94|l is 

d r {^gg u g rr SF tr + ^gg^g rr SF^ r ) + d z W=gg t *g" > 5F^)]F^ ) + 
d r {V=gg u g rr SF^ r + ^g^g rr 6F tr ) + d z {^gg^g zz 5F^ z ) + 
dt (V^gg 4 " p g tt SF^ t + V=^<V*«fy) 

where the terms involving SF tz have canceled. The terms involving 5F tr may be combined to give 



F (0) - 

C 4>r — U - 



-gg^g^SFr — 



d f-9**\ „(0) 



dr V 9 tt ) r0 



f: 



(Cl) 



(C2) 
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The terms involving SF^r combine to give 

f=g 8 



Or 



9rr9<t»t> dr 



F 



(0) 



d 



r<p 



g**\ d 
dr \ g u ) dr 



Here det(i, <fi) = gu9c/></> — 9td>- Lastly, the terms involving 6F$t and SF^ yield 



-g I d 2 Sr g zz d 2 Sr 



(0)\2 



det(i, (f>) V dt 2 g u dz 2 
To first order in the coordinate perturbation, the perturbed Maxwell fields are 

d_ 

dr 



SF r . 



(C3) 



(C4) 



(C5) 



SF rt = ^f(r) " 



dt 



dr 



9 tt (r) r * ( \ 



(C6) 



x F _ 9<5r p(0) M . 
* ~ ~dt ^ 



and 



SF, 



dz r <t> (ry 



Substituting these expressions and collecting terms, we arrive at 



1 



det(t, c 



d 2 Sr g zz d 2 Sr 



dt 2 g u dz 2 



-gdr 



-g d 



T^F^) 

9^9rr dr v 



F 



(0) 



r<p 



9 U 9 rr 



d 



.-10 



dr \ g ft 

If we multiply by det(^>, t)g tt we obtain equation (|10D|> . 



dt dr \ g u 



= 0. 



(C7) 



(C8) 



(C9) 



APPENDIX D 



This www.arxiv.org manuscript incorporates several minor corrections that will appear as an erratum to the paper 
as published in the December 2004 issue of Physical Review D. 

1) The second term, RHS, Eq. (100) [and the analogous term in Eq. (C9)] was spurious and has been deleted. 

2) An algebraic error was corrected in the RHS of Eqns. (181) and (182): the normalization increased by a factor 
2; and r _4 [r 4 ] — * r _1 [r] in front of [inside] the integral. The intermediate equation (179) also received a correction. 

3) In Section V.B, a final step has been added to the derivation of the fast mode amplitude [Eq. (166b)]. The claim 
[based on Eq. (157)] that 6R max can be set to zero was incorrect: the integral for the background field energy in Eq. 
(157) in fact converges at large radius. As a result, the LHS of Eq. (166) gains a term — SR max / R msix , and the final 
step Eq. (166b) is required. 

In addition there were a handful of (non-propagating) corrections to signs and indices. 
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